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ABSTRACT 



Through the Optimum Path Aircraft Routing System (OPARS), 
Fleet Numerical Oceanography Center (FNOC) has committed it- 
self to providing a computerized flight planning service 
remotely accessible via dial-up communications lines. The 
question arises as to whether the proposed number of tele- 
phone lines will be adequate to provide a level of service 
previously provided by the Lockheed Jetplan system. This 
study provides a detailed analysis of the response delays for 
the OPARS flight plan system. In addition, estimates are 
given of communication requirements when various levels of 
demand prevail, and under conditions in which the FNOC computer 
is busy or idle. 
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NOTATIONS AND ABBREVIATIONS 



A/B/m m-Server queueing system where A represents the 

the interarrival distribution and B represents the 
service distribution 

CM Central Memory 

CPU Central Processor Unit 

E Denotes r stage Erlang distribution 

E[T(Setup)] Expected service time of IRG sub-system 

E[T(Input queue)] Expected service time of Input queue sub-system 

E[T(Run)] Expected service time of execute sub-system 

E[T(Service) ] Expecter service time for OPARS system 

PNOC Fleet Numerical Oceanography Center 

G Denotes general distribution 

X Q Arrival rate of 0PAR3 jobs 

Arrival rate of higher priority jobs 

y Service rate 

M Denotes Exponential distribution 

N(t) Number of customers in the system at time t 

N Long run average number of customers in the system 

OPARS Optimum Path Aircraft Routing System 

p. Erlang formula state probabilities for states 

i = 1,2, ... ,11 

p utilization factor 
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I. SUMMARY 



In the proposed OPARS flight plan system, users at 
37 remote terminals (Figure 1) connect directly to the 
Fleet Numerical Oceanography Center (FNOC) computer system 
via 11 telephone lines. The total amount of time that a line 
is busy, the total OPARS service time, can be modeled as 
the sum of three sub-system service times: 

1. An interactive program setup time T(setup); 

2. A delay that the OPARS program experiences in 
the FNOC computer's input queue T( input queue); 

3. An OPARS program run time T(run). 

The expected values of these individual sub-system service 
times can be computed and then summed to provide a mean 
or expected service time for the entire OPARS request as 
follows : 

E[T( service) ] = E[T(setup)] + E[T(input queue)] + E[T(run)]. 
From this. Erlang's formula equation (l) can be employed to 
compute individual long-run state probabilities (the 
probability that i lines are busy), p^, and ultimately the 
probability, P^j of having all eleven lines busy. 

In support of the analytic model, a computer program 
was written to simulate the entire OPARS flight plan request 
process, from initial dial-up through program completion. 

The simulation also includes the arrival and servicing of other 
FNOC computer programs which interfere with the processing 
of the OPARS program. 
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Figure 1. Remote Terminal Sites 
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Both the analytic model and the simulation, computed 
state probabilities for OPARS demand rates of 2, 4, 6, ...20 
per hour and under conditions in which the FNOC computer is 
busy or idle. The results clearly indicate that the 11 
telephone lines can handle demand rates up to three times 
as great as those currently being experienced by the Lockheed 
Jetplan system before p , probability of all lines 

busy, exceeds, .05- 
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II. INTRODUCTION 



Fleet Numerical Oceanography Center (FNOC) is currently 
testing and evaluating a computerized flight plan system, 
referred to, for short, as OPARS. This sytem, developed to 
replace the Lockheed Jetplan flight plan sytem, provides users 
at remote sites with direct access to the FNOC computer 
via 11 telephone lines. The purpose of this study is to 
determine if the intended number of telephone lines would 
be adequate to ensure a low probability of having all lines 
busy. 

The number of lines busy at any time t, {N( t) | t>o} can be 
modeled as a Birth-Death process with inter-arrival times 
that are assumed to be independent and exponentially 
distributed but with service times that clearly are not. 
Fortunately, in a communications system characterized by 
calls that arrive in a stationary or time-homogenous Poisson 
Process of rate X and vie for n lines with no queue (blocked 
calls lost), the long run probability of having i lines busy, 
Pj_, can be computed from Erlang’s formula. 



r _ (At) /i! 

p i ~ i H Z 

2 (Xt)/j ! 

j=0 



0<i<N 



i>N 



( 1 ) 



0 



This formula has the surprising feature of being valid for 
any service time distribution. Indeed, given the above 
assumption, one need only obtain the mean service time, t, 



11 



in order to calculate p. . The majority of the effort in 
this study was to characterize the OPARS system in such a 
way so that a mean service time could be derived under various 
operating conditions. 

This paper begins with some introductory material 
concerning the PNOC computer systems, the OPARS flight plan 
system and historical usage of the Lockheed Jetplan system. 
Section IV is a detailed description of the analytic model 
used and a discussion of the supporting siumulation program. 
Finally, Section V contains the tabled results and an analysis. 

Notation used in this paper is consistent with that found 
in Queueing Theory literature and will be explained whenever 
introduced. In addition, a listing and discription of 
all notation and abbreviations can be found on page seven. 
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III. BACKGROUND 



A. FNOC COMPUTER SYSTEM 

FNOC currently has four computers in service at their 
Monterey facility, but only one, referred to as HAL is 
accessible to the remote OPARS terminals. HAL is a CDC-6600 
computer with two central processors and 330 K of octal 
central memory. HAL is accessed within FNOC by approximately 
twenty-five interactive terminals as well as by a batch 
job system that routinely processes hundreds of programs a 
day. Since HAL alone is involved in providing service to 
OPARS remote terminals, we will only consider it in the 
discussion. 

1 . HAL Batch System 

HAL's batch system is composed of a priority ranked 
"first-in, first-out" input queue, which can be considered 
exterior to the computer, and a priority ranked execute 
queue from which jobs are selected by the scheduler for 
processing (Figure 2). When a job enters the system, it 
first enters the input queue, at which the job's user-assigned 
priority establishes its initial position. While in the 
queue the job slowly accrues additional "wait time" 
priority which ensures that even the lowest priority jobs 
are eventually run. The largest number of daily jobs are 
restricted to priorities of 1, 2 or 3 (3 being the highest) 
but priorities up to 77 can be assigned. The priority of 
OPARS jobs is fixed automatically at 60 which immediately puts 
it ahead of all but a few other jobs in the queue. 
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When space is available within the computer the scheduler 
removes the first job (with respect to priority then to 
time-of-arrival ) from the input queue and moves it to the 
execute queue. Upon entering the execute queue the job 
loses all of its previously gained "wait-time" priority but 
immediately begins to gain it back as the job waits to be 
processed. In this queue the jobs with the higher assigned 
priorities gain this additional priority faster than those 
with low assigned priority. This mechanism results in higher 
priority jobs getting selected for processing sooner than 
low priority jobs. When selected by the scheduling routine, 
the program is moved into central memory and it is processed 
for a unit length of time or until it attempts to access 
a device (disk, tape, etc.) which is unavailable. When this 
event occurs the job is removed form central memroy and 
returned to the execute queue, having its accured priority 
reset to zero. The process continues in this manner until 
the program is completed and it is transferred to the out- 
put queue. 

2 . Current FNOG Workload 

From available central processor and central memory 
utilization profiles (Figure 3) it can be readily seen that 
resource demands on HAL can be separated into two states: 

A Busy or High Demand state in which the computer is operating 
at maximum capacity, and an Idle or Low Demand state in which 
most of the resources are immediately available. The High 
Demand state is characterized by a full execute queue, and a 
lengthy input queue, while the Low Demand state exhibits an 
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Figure 2. FNOC Computer System 





execute queue with only an occasional job in it and an in- 
put queue which is empty. Fortunately the transition 
period between these states is quite brief and therefore 
a transition state is not needed. 

The High Demand period coincides, not surprisingly > 
with normal working hours except that it extends well into 
the evening, often as late as midnight. Any jobs entering 
the system during this period must be delayed for a time 
in the input queue before being run. During the idle period, 
on the other hand, jobs pause only momentarily in the input 
queue before being run. 

B. OPARS FLIGHT PLAN SYSTEM 

The OPARS flight plan system provides the user with 
an optimum (with respect to time) route of flight given forecast 

weather and user inputs such as aircraft type and, take off 
weight. For flights within areas that have a defined 
routing structure the optimization is fairly easy because 
only a few routes are candidates for the optimum. Over 
water or point to point flight plans are more difficult 
because there is, in theory, an uncountably infinite number 
of available routes. 

There are two major sections of the OPARS program: an 

interactive portion and the optimization program. When a 
user at a remote site requests an OPARS flight plan one 
first sets up the program with the Interactive Response 
Generator (IRG). Using a query and response technique, the 
IRG prepares the main program by obtaining required 
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i 2 



Figure 3. CPU and CM Utilization Profiles 



information from the user. During this session the IRG 
is not checking input information for validity, but only 
for format. For example, an entry of ABCE, as the four- 
letter identification code for the destination airfield, 
would be accepted even though no such airfield code 
exists; When all of the required information is entered, 
the IRG allows the user to review and change any input 
if necessary. Then, upon command, the program is put into 
the HAL computer batch system. 

As discussed in the previous section, the OPARS job 
enters the input queue along with all other FNOC jobs and 
awaits its turn. Hereafter the OPARS program is handled 
as any other batch program with one exception. Upon 
completion the flight plan is returned automatically to the 
terminal where the request was entered. If for some 
reason the line has been disconnected^ the flight plan is 
placed into an output file from which it can be retrieved 
later. 

C. HISTORICAL USAGE/PROJECTED DEMAND 

In an effort to evaluate expected demand for the OPARS 
flight plan, several methods were employed. They were: 

1) A linear regression model of historical usage by 
the Navy to estimate demand through 1981; 

2) expectations of remote site users; 

3) Lockheed's records of previous use. 

None of these methods alone provided the required level of 
accuracy and, only by combining the three, could a reasonable 
estimate be made. 
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The Navy's own records provided the best overall 
information. This information, shown in Figure consists 
of monthly totals of Navy requests for the Jetplan. The 
trend is clearly an increase in usage and a simple linear 
regression yielded the projected demands shown in Figure 5 . 
These monthly totals, unfortunately did not provide any 
information as to how requests varied during the day. 
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Figure 5. Linear Regression of Usage Data 
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In order to obtain the distribution of daily requests, 
expected demand information was requested from all remote 
site users. The resulting information suffered from 
imprecision and a lack of completeness (only about half of 
the remote sites responded) . 

As a final effort, an attempt was made to secure actual 
daily records from the Lockheed Corporation. The Lockheed 
personnel were unable to provide me with actual data, but 
over the course of several interviews a picture of the 
daily demand distribution was pieced together. 

The results of these data collection efforts appear in 
Figure 7. These graphs show the projected demand for 
July 1980 and July 1981 distributed using the composite 
results of Lockheed information and user expectations. 



OPARS 

REQUESTS 
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Local Hours 

Figure 6. Predicted Demand Rates for OPARS 



IV. THE MODEL 



We start the modeling discussion by reviewing 
some basic assumptions about the system. 

First and foremost is the assumption concerning 
time-homogeneous Poisson arrivals for OPARS requests. 
Without this assumption, the Erlang formula (1) would not 
be valid and a much more complex model would be required. 

Secondly, the service times for the three sub-systems 
are assumed to be mutually independent. 

Finally, as discussed in the last section, the 
conditions for the computers workload will be divided in- 
to two distinct periods. These are the High Demand period 
with a nonempty input queue, and the Low Demand period 
where the OPARS requests bypass the input queue and go 
directly into the computer. 

A. THE ANALYTIC MODEL 

Since the main thrust of the analytic effort is to 
obtain an overall mean service time, an immediate 
simplification would be to somehow break the OPARS system 
into two or more sub-systems whose service times would 
be more easily obtained. With this in mind, we start by 
separating the service into the interactive sub-system 
and the batch sub-system. The latter is still rather 
formidable because of the nature of the input and 
execute queues and the complexity of the scheduling 
routine;. The final step is to reduce the batch system 



into an input queue sub-system and an execution sub-system. 
The reasons for this may not be readily apparent as there 
are other modeling techniques such as an M|G|m multiserver 
queueing system with priorities, which would handle the 
system as a whole. Such a system would require that each 
job priority have a mean service time. Unfortunately, there 
is no correlation between a jobs priority and its core 
and central processor requirements. Additional unrealistic 
assumptions would be required which would, in the end, 
reduce the validity of the result. 

1. T he Interactive Sub-System 

First to be considered is the interactive sub-system. 
An estimate for the mean time for this sub-system is 
readily comput ble from the available, albeit scarce, data. 
This data(Figure 7) was obtained by timing FNOC technicians 
on trial OPARS runs. The mean of these data is 5-1 minutes 
and will be assumed to be the same during both High and 
Low Demand periods. 

2 . Input Queue Sub-System 

The next sub-system to model is the input queue sub- 
system. Under low demand conditions the mean service time 
is, by definition, equal to zero. During the High Demand 
period it is not quite so simple. Again, by previous 
definition, the High Demand period has a continuously non- 
empty input queue. This very simple assumption suggests that 
the input queue may be modeled as a single-server queueing 
system by itself, the server being the first position in the 
queue. One must merely postulate the interarrival process 
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and the service process, the latter representing the times 
between successive departures from the input queue. If we 

x 

can now model the arrival rate as Poisson and the service 
distribution as exponential we have an M/M/1 queueing 
system ranked by priority for which there is a simple 
closed form solution for expected waiting times. The 
interested reader is urged to consult Griffin [2] for the 
derivation of this result. 

There was no data available to determine the precise 
manner in which OPARS jobs enter the input queue and so a 
Poisson arrival rate is as plausible as any other. However 
the service distribution (the interdeparture times from the 
input queue) can be measured by observing a real time display 
of HAL's input queue and timing departures. The resulting 
data (Figure 8) and histogram (Figure 9) are surprisingly 
close to the exponential distribution. In fact, the Chi- 
squared goodness of fit test resulted in a test statistic of 
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Figure 8. Input Queue Interdeparture Time 
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Figure 9. Histogram of Input Queue Interdeparture Time Data 
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0. 533 j which indicates that the data is very close to being 
exponential with Aq = 53-52 per hour. 

The final consideration is the priority ranking. 

As stated earlier an OPARS job is assigned a priority of 6 0 
while FNOC jobs can have priorities from 1 to 77- However, 
the batch jobs of lower priority can be ignored and we are 
left with only higher priority jobs to consider. Available 
data (Figure 10) showed higher priority jobs arriving at the 
rate of 3-^5 per hour. The one remaining assumption is that 
these higher priority jobs also arrive according to a 
Poisson Process. 

The expression for OPARS expected input queue is 

then, 

E[T( INPUT QUEUE)] = (2) 

y ( 1 - fJn ( i - A nj 
U U 

where u is the service rate, y is the OPARS arrival rate and 
A is the higher priority job arrival rate. This expression 
is valid only when: 

A h + X y > 1. 

U 

This input queue waiting time was computed for OPARS arrival 
rates of 2, 4, 6,..., 20 per hour with results listed in 
column three of Figure 16. 

3. The Execution Sub-Systems 

In the third and final section, the execution 
sub-system, we will characterize the OPARS program run time as 
well as delays due to other factors. 
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Because of the central memory requirements of the OPARS 
program (150K) only two OPARS jobs are allowed in the computer 
at one time. Any additional requests arriving during such 
a state would be required to wait in the input queue and would 
be passed over by the scheduling routine as it selected jobs 
for processing. Because the two OPARS jobs can be processed 
simultaneously, a queueing system of the form GI/G/2 is 
suggested. This type of system is uncomfortable to work 
with, so simplifying assumptions will be made. 

First, because of a lack of data to the contrary, 
we will again ssume Poisson arrivals of OPARS jobs. It 
now remains to describe the service distribution. Unfortunately, 
even with Poisson arrivals, multi-server systems with any but 
exponential service times, are difficult to analyze math- 
ematically . 

Figures 11 and 12 contain actual OPARS run-time data 
for Low Demand and High Demand periods respectively. These 
data reflect the time an OPARS job spends in the system 
from it’s transfer to the execute queue until it’s completion. 
Histograms of these data are in Figures 13 and 14. and are 
clearly not exponential . However, if the data can be 
considered to be from an Erlang distribution, then the nomo- 
graph in Figure 15 (see Hillier and Leiberman [4]) can be 
used to obtain an approximate value for N, the mean number 
of OPARS jobs in the system. From this the mean service time 
can be computed via Little's result: 

E[T(run) ] = - 
^ 0 

where A 0 is, as before, the arrival rate for OPARS jobs. 
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Figure 13. Histogram of OPARS Program Run time Data; Low Demand 
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Figure 15. Nomograph of the Long Run Average 
Number of Customers, N, in an M/E^/2 Queue 
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Figure 17. Idle System Summary (Time in Minutes) 



The run data was parameterized as an Erlang 



distribution with the following results: 
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Where 




is 



the mean run time. 



Finally, the expected execution service time, 
E[T(runs)], was computed for a series of OPARS arrivals, 
as before, with the results in column 4 of Figures 16 . and 1'7. 

4 . Summar y 

In summary we have now succeeded in representing 
the total OPARS service time as the sum of three, readily 
computable sub-system times (Figure 19) : 

E[T( service) ] = E[T(setup)] + E[T(input queue)] +E[T(run)]. (4) 

Erlang's formula can now be employed to compute expected 
state probabilities for lines busy. 
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Figure 19- OFARS System Diagram 



B . THE SIMULATION 



As an alternative solution technique, two computer 
programs were written to simulate the OPARS system, one each 
for High and Low Demand computer states. Output of the 
simulations were, as in the analytic' model, state probabilities 
for lines busy under various demand states. 

There are two key differences in the logic flow for the 
two simulations. First OPARS requests arriving during the 
busy state are sent to the input queue where they must wait 
to be served. Under idle conditions jobs go directly to 
the computer. Secondly, in the busy state, higher priority 
jobs are being created which interfere with the processing 
of OPARS jobs. The following is a list of the various 
distributions and their parametric values used in the 
simulat ion. 

OPARS Interarrival time: EXP (\= 2 , 4 , 6 , . . . 20/HR) 

High Priority Jobs Interarrival time:.. EXP(x = 3 . 45 /HR) 

IRG service times: NORMAL (y= 5-1 Min., 0= 1) 

Input queue Interdeparture times : EXP(x= 53. 5/HR) 

OPARS run times » Idle: GAMMA (X= .721, r = 4.34) 

OPARS run times.. Busy: GAMMA (A= .701, r = 3.1) 

Figure 20: Probability Distribution Used in Simulation 

Program 

The simulations were run for 48 hours for each OPARS 
demand rate and under each computer state. 
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The only significant factor included in the simulations 
but not in the analytic model had to do with a retrial 
population. In the event of all lines busy, the analytic 
model assumes arriving requests disappear, but in reality, 
people, when faced with a busy signal, hang up and try 
again. The simulation retries all balked attempts after 
ten minutes. 
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V. RESULTS AND ANALYSIS 



Erlang formula (1) state probabilities, pi, were 
computed for OPARS request rates of 2, 4,6,..., 20 per hour 
under High and Low Demand conditions. Appendix A contains 
the FORTRAN program used to calculate these state probabilities 
given demand rate and mean service time, E[T( service )] . 
Similarly Appendices D and C contain the SIMSCRIPT programs 
which simulate the OPARS system. 

The results of both methods are in Figures 2. and 22. 

The state probabilities are suprisingly close (differing 
by less than 2 % in most cases) and clearly show that with 
eleven operating telephone lines, demand rates must be far 
in excess of projected requirements (Section III.C), before 
the probability of all lines being busy is of concern. 

One should be aware, however, that in spite of the 
closeness of the results these two solutions merely serve 
to verify each other. Because the OPARS system is not fully 
functional at the time of this report, there is no way 
to validate either method against the actual system. 

In addition, if any significant changes are made to 
either the OPARS system or to the FNOC computer system, 
a new evaluation would be required. 
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Figure 21. Computed Erlang Formula State Probabilities 



STATE 2/HR 4/HR 6/HR 8/HR 10/HR 12/HR 14/HR 16/HR 18/HR 20/HR 

0 .6646 .4238 .2894 .1890 .1313 .0815 .0243 .0175 .0069 .0006 
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Figure. 22 . Simulation Program State Probabilities 



APPENDIX A 



COMPUTER PROGRAM TO CALCULATE ERLANG FORMULA STATE PROBABILITIES 



DIMENSION PROB (12,10), Serve -(10) 

ZL = .0001 
NN = 1 

5 READ (5,10) ( SERVE ( I ) , I = 1,10) 

10 FORMAT (IF10. 5) 

LINE = 12 

DO 60 J = 1,10 

RHO = 2.0 * J/SERVE( J ) 

DENOM = 0.0 
DO 30 K = 1 , LINE 
FACT = 1.0 
DO 20 L = 1 , K 
FACT = FACT * (L-l) 

IF (L . EQ. 1) FACT = 1.0 
20 CONTINUE 

DO 50 I = 1 , LINE 

FACTNU =1.0 

DO 40 L = 1.1 

FACTNU = FACTNU * (L-l) 

IF(L .EG. 1) FACTNU = 1.0 
40 CONTINUE 

PROB (I,J) = ( RHO* * ( 1-1 ) /FACTNU ) /DENOM 
IF ( PROB ( I , J ) .LT. ZL) PROB(I,J) =0.0 
50 CONTINUE 
60 CONTINUE 

WRITE (6,90) 

DO 70 I = 1, LINE 
M = 1-1 

WRITE( 6,65) M, ( PROB ( I , J ) , J = 1,10) 

65 FORMAT ( 2X, 12 , 10F11 . 6 ) 

70 CONTINUE 

IF (NN.EQ.2) STOP 
NN = 2 
GO TO 5 

90 FORMAT Cl) 

END 
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APPENDIX B 



SIMULATION PROGRAM: HIGH DEMAND STATE 
PREAMBLE 

NORMALLY MODE IS INTEGER 

EVENT NOTICES INCLUDE HIGH . PRI .JOB , OPARS.JOB, RETRY, 
INPUT. QUEUE, 

TEMPORARY ENTITIES 

EVERY JOB HAS A TIME . OF . ENTRY AND A PRIORITY AND MAY 
BELONG TO THE QUEUE 

DEFINE QUEUE AS A FIFO SET RANKED BY PRIORITY 
THE SYSTEM OWNS THE QUEUE 

DEFINE IRG, THRUPUT, QUEUE. TIME AND TIME. OF ENTRY AS 
REAL VARIABLES 

DEFINE NO. REQUESTS, NO. ATTEMPTS, PRIORITY, LINES. 

BUSY, RUN, NO. BATCH, SYSTEM, BUSY AND EMPTY AS VARIABLES 
DEFINE SUMMARY AS A 2-DIMENSIONAL, REAL ARRAY 
ACCUMULATE STATE. PROB(0 to 11 BY 1) AS THE HISTOGRAM OF 
LINES, BUSY 

ACCUMULATE QUEUE . STATS ( 0 TO 11 BY 1) AS THE HISTOGRAM 
OF N. QUEUE 

TALLY AVE. QUEUE. TIME AS THE AVERAGE OF QUEUE. TIME 



MAIN 

RESERVE SUMMARY ( * , * ) AS 12 BY 10 
LET RUN = 1 
LET BUSY = 1 

SCHEDULE AN OPARS.JOB NOW 

SCHEDULE A STATISTICS IN 2830 MINUTES 

START SIMULATION 

END 

EVENT OPARS.JOB 

SCHEDULE AN OPARS.JOB IN EXPONENTIAL . F ( 30 . /RUN , 4 ) MINUTES 

ADD 1 TO NO. ATTEMPTS 

IF LINES. BUSY EQUALS 11 

SCHEDULE A RETRY IN 10 MINUTES 

RETURN 

OTHERWISE 

ADD 1 TO LINES. BUSY 

ADD 1 TO NO. REQUESTS 

LET IRG = NORMAL. F(5. 1,1. 0,5) 

SCHEDULE AN INPUT. QUEUE IN IRG MINUTES 
RETURN 

END 

EVENT HIGH. PRI. JCB 

SCHEDULE A HIGH. PRI. JOB IN EXPONENTIAL . F( 17 • ^ , 3) MINUTES 

CREATE A JOB 

LET PRIORITY (JOB) = 2 

FILE JOB IN QUEUE 

RETURN 

END 
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EVENT RETRY 

IP LINES. BUSY EQUALS 11 

ADD 1 to NO. ATTEMPTS 

SCHEDULE A RETRY IN 10 MINUTES 

RETURN 

OTHERWISE 

ADD 1 TO NO. REQUEST 
ADD 1 TO LINES. BUSY 
LET IRG = NORMA.F(5.1,1.0,5) 

SCHEDULE AN INPUT. QUEUE INIRG MINUTES 
RETURN 

END 

EVENT INPUT. QUEUE 

IP SYSTEM EQUALS EMPTY AND NO. BATCH IS LESS THAN 2 
SCHEDULE A JOB . COMPLETION IN GAMMA . F ( 4 . 4 3 , 3 • 1 , 7) MINUTES 
ADD 1 TO NO. BATCH 
OTHERWISE 
CREATE A JOB 

LET TIME. OP. ENTRY (JOB) = TIME.V 

LET PRIORITY (JOB) = 1 

FILE JOB IN QUEUE 

REGARDLESS 

RETURN 

END 

EVENT EXECUTION 

IP SYSTEM EQUALS BUSY 

SCHEDULE AN EXECUTUION IN EXPONENTIAL . F( 1 . 12 , 2) MINUTES 

REGARDLESS 

IF N. QUEUE EQUALS 0 

RETURN 

OTHERWISE 

IF PRIORITY (F. QUEUE) = 2 

REMOVE FIRST JOB FROM QUEUE 

DESTROY THE JOB 

RETURN 

OTHERWISE 

IF NO .BATCH . EQUALS 2 

RETURN 

OTHERWISE 

REMOVE FIRST JOB FROM QUEUE 

LET QUEUE. TIME= (TIME.V -TIME . OF. ENTRY ( JOB) * 1440. 
DESTROY THE JOB 
IF SYSTEM EQUALS BUSY 

SCHEDULE A JOB . COMPLETION IN GAMMA . F( 6 . 06 , 4 . 34 , 6 ) MINUTES 
OTHERWISE 

SCHEDULE A JOB . COMPLETION IN GAMMA . F( 4 . *13 , 3 . 1 , 7) MINUTES 

REGARDLESS 

ADD 1 TO NO. BATCH 

RETURN 

END 
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EVENT JOB. COMPLETION 

SUBTRACT 1 FROM LINES. BUSY 
SUBTRACT 1 FROM NO. BATCH 
IF N. QUEUE EQUALS 0 
RETURN 
OTHERWISE 

IF SYSTEM EQUALS BUSY 
CANCEL THE EXECUTION 
REGARDLESS 

SCHEDULE AN EXECUTION NOW 
RETURN 

END 

EVENT STATISTICS 

START NEW PAGE 

PRINT 1 LINE WITH RUN*2 THUS 
ARRIVAL RATE = **HOUR 

SKIP 3 OUTPUT LINES 

PRINT 1 LINE THUS 

PR( BUSY LINES) PR(INPUT QUEUE) 

SKIP 2 OUTPUT LINES 
FOR I = 1 TO 12, DO 

PRINT 1 LINE WITH 1-1, STATE . PROS ( I ) /2 AND QUEUE. 

STATS ( I) /2 THUS 
** ***** ***** 

LOOP 

SKIP 5 OUTPUT LINES 

PRINT 3 LINES WITH AVE. QUEUE. TIME, NO. ATTEMPTS AND NO. 
REQUESTS THUS 

AVERAGE OPARS QUEUE TIME = **.**** 

NO . OPARS ATTEMPTED **** 

NO. OPARS COMPLETE **** 

FOR I = 1 to 12, DO 

LET SUMMARY (I, RUN) = STATE . PROS ( I ) /2 
LOOP 

SCHEDULE A STATISTICS IN 2880 MINUTES 

RESET TOTALS OF LINES. BUSY, N. QUEUE AND QUEUE. TIME 

IF RUN EQUALS 10 

GO TO FINAL. STATS 

OTHERWISE 

ADD 1 TO RUN 

LET NO. ATTEMPTS ' = 0 

LET NO. REQUESTS = 0 

RETURN 

’ FINAL. STATS r 
START NEW PAGE 
PRINT 5 LINES THUS 
EMPTY SYSTEM SUMMARY 

STATE 2/HR 4/HR 6/HR 8/HR 10/HR 12/HR 14/HR 16/HR 18/HR 20AIR 
FOR I = 1 to 12, DO 

PRINT 1 LINE WITH 1-1, SUMMARY (1,1), SUMMARY (1,2), 
SUMMARY ( I , 3 ) , SUMMARY (1,4), SUMMARY (I , 5 ) , SUMMARY (1,6) , 
SUMMARY ( I , 7 ) j SUMMARY (1,8) , SUMMARY (1,9) , SUMMARY (1,10) 

** **** **** **** **** **** **** **** **** **** ' ** 

• ••••••••• 

LOOP 

STOP 

END 46 






APPENDIX C 



SIMULATION PROGRAM: LOW DEMAND 
PREAMBLE 

NORMALLY MODE IS INTEGER 

EVENT NOTICES INCLUDE HIGH . PRT . JOB , OP ARS .JOB , RETRY , INPUT . 
QUEUE, 

TEMPORARY ENTITIES 

EVERY JOB HAS A TIME . OF . ENTRY AND A PRIORITY AND MAY 
BELONG TO THE QUEUE 

DEFINE QUEUE AS A FIFO SET RANKED BY PRIORITY 
THE SYSTEM OWNS THE QUEUE 

DEFINE IRG,THRUPUT, QUEUE. TIMS AND TIME . OF . ENTRY AS 
REAL VARIABLES 

DEFINE NO. REQUESTS; NO. ATTEMPTS , PRIORITY, LINES. BUSY, 

RUN, NO. BATCH, SYSTEM, BUSY AND EMPTY AS VARIABLES 
DEFINE SUMMARY AS A 2-DIMENSIONAL, REAL ARRAY 
ACCUMULATE STATE. PROB(0 TO 11 BY 1) AS THE HISTOGRAM 
OF LINES. BUSY 

ACCUMULATE QUEUE . STATS .( 0 TO 11 BY 1) AS THE HISTOGRAM 
OF N. QUEUE TALLY AVE . QUEUE .TIME AS THE AVERAGE OF QUEUE. 
TIME 

END 

MAIN 

RESERVE SUMMARY ( * , *) AS 12 BY 10 
LET RUN = 1 
LET SYSTEM = 1 
LET BUSY = 1 

SCHEDULE AN OP ARS . JOB NOW 

SCHEDULE AN EXECUTION IN EXPONENTIAL. F( 1, 12, 2) MINUTES 

SCHEDULE A HIGH. PRI. JOB IN EXPONENTIAL . F( 17 . 4 , 3) MINUTES 

CREATE A JOB 

LET PRIORITY (JOB) = 2 

FILE JOB IN QUEUE 

SCHEDULE A STATISTICS IN 2880 MINUTES 
START SIMULATION 

END 

EVENT OP ARS . J OB 

SCHEDULE AND OPARS.JOB IN EXPONENTIAL . F( 30 . /RUN, 4 ) MINUTES 

ADD 1 TO NO. ATTEMPTS 

IF LINES. BUSY EQUALS 11 

SCHEDULE A RETRY IN 10 MINUTES 

RETURN 

OTHERWISE 

ADD 1 TO LINES. BUSY 

ADD 1 TO NO. REQUESTS 

LET IRG = NORMA. F(5. 1,1. 0,5) 

SCHEDULE AN INPUT. QUEUE IN IRG MINUTES 
RETURN 

END 

EVENT HIGH. PRI. JOB 

SCHEDULE A HIG.PRI.JOB IN EXPONENTIAL. F( 17 • 4 , 3) MINUTES 
CREATE A JOB 
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APPENDIX C 



SIMULATION PROGRAM: LOW DEMAND STATE CONT'D 

LET PRIORITY( JOB) = 2 
PILE JOB IN QUEUE 
RETURN 

END 

EVENT RETRY 

IF LINES. BUSY EQUALS 11 

ADD 1 TO NO .ATTEMPTS 

SCHEDULE A RETRY IN 10 MINUTES 

RETURN 

OTHERWISE 

ADD 1 TO NO. REQUESTS 

ADD 1 TO LINES .BUSY 

LET IRG = NORMAL. F(5. 1,1. 0,5) 

SCHEDULE AN INPUT. QUEUE IN IRG MINUTES 
RETURN 

END 

EVENT INPUT. QUEUE 

IF SYSTEM EQUALS EMPTY AND NO. BATCH IS LESS THAN 2 
SCHEDULE A JOB . COMPLETION IN GAMMA .F(4. 43,3-1,7) MINUTES 
ADD 1 TO NO .BATCH 
OTHERWISE 
CREATE A JOB 

LET TIME. OF. ENTRY (JOB) = TIME.V 

LET PRIORITY (JOB) = 1 

FILE JOB IN QUEUE 

REGARDLESS 

RETURN 

END 

EVENT EXECUTION 

IF SYSTEM EQUALS BUSY 

SCHEDULE AN EXECUTION IN EXPONENTIAL. F ( 1. 12,2) MINUTES 

REGARDLESS 

IF N. QUEUE EQUALS C 

RETURN 

OTHERWISE 

IF PRIORITY (F. QUEUE) = 2 

REMOVE FIRST JOB FROM QUEUE 

DESTORY THE JOB 

RETURN 

OTHERWISE 

IF NO. BATCH EQUALS 2 

RETURN 

OTHERWISE 

REMOVE FIRST JOB FROM QUEUE 

LET QUEUE. TIME = (TIME.V - TIME . OF . ENTRY (JOB ) ) * 1440 . 
DESTORY THE JOB 
IF SYSTEM EQUALS BUSY 

SCHEDULE A JOB . COMPLETION IN GAMMA . F( 6 . 06 , 4 . 34 , 6 ) MINUTES 
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SIMULATION PROGRAM: LOW DEMAND STATE CONT'D 
OTHERWISE 

SCHEDULE A JOB . COMPLETION IN GAMMA. F( 4 . i| 3 , 3 . 1 , 7 ) MINUTES 

REGARDLESS 

ADD 1 TO NO. BATCH - 

RETURN 

END 

EVENT JOB. COMPLETION 

SUBTRACT 1 PROM LINES. BUSY 
SUBTRACT I FROM NO. BATCH 
IF N. QUEUE EQUALS 0 
RETURN 
OTHERWISE 

IF SYSTEM EQUALS BUSY 
CANCEL THE EXECUTION 
REGARDLESS 

SCHEDULE AND EXECUTION NOW 
RETURN 

END 

EVENT STATISTICS 

START NEW PAGE 

PRINT 1 LINE WITH RUN *2 THUS 
ARRIVAL RATE = **/HOUR 

SKIP 3 OUTPUT LINES 

PRINT 1 LINE THUS 

PR( BUSY LINES) PR( INPUT LINES) 

SKIP 2 OUTPUT LINES 
FOR I = 1 TO 12, DO 

PRINT 1 LINE WITH 1-1, STATE. PROBLE(I) . 2 AND QUEUE. STATS 
(I)/2 THUS 

## .#**** .#***# 

LOOP 

SKIP 5 OUTPUT LINES 

PRINT 3 LINES WITH AVE . QUEUE . TIME .NO . ATTEMPTS AND NO. 
REQUESTS THUS 

AVERAGEOPARS QUEUE TIME = **.**** 
no.opars ATTEMPTED *##* 

NO.OPARS completed **** 

FOR I = 1 TO 12, DO 

LET SUMMARY (I, RUN) = STATE . PROB( I ) /2 
LOOP 

SCHEDULE A STATISTICS IN 2880 MINUTES 

RESET TOTALS OF LINES .BUSY, N. QUEUE AND QUEUE. TIME 

IF RUN EQUALS 10 

GO TO FINAL. STATS 

OTHERWISE 

ADD 1 TO RUN 

LET NO. ATTEMPTS = 0 

LET NO .REQUESTS = 0 

RETURN 

'FINAL. STATS’ 

START NEW PAGE 
PRINT 5 LINES THUS 
BUSY SYSTEM SUMMARY 
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STATE 2/HR 4/HR 6/HR 8/HR 10/HR 12/HR 14/HR 16/HR 18/HR 20/HR 
FOR I = 1 TO 12, DC 

PRINT 1 LINE WITH 1-1, SUMMARY (1,1), SUMMARY (1,2), 
SUMMARY (1,3) ,SUMMARY( I , 4) ,SUMMARY(I,5) , SUMMARY (1,6) , 
SUMMARY (1,7), SUMMARY (1,8), SUMMARY (1*9), SUMMARY 
(1,10) THUS 

## #### #### #### **** #### #### #### #### ###* 

• ••••••••• 

LOOP 

STOP 

END 
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